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Abstract. A system of two biased, mutually exclusive random walkers on an infinite 
i-G . ID lattice is studied whereby the intrinsic bias of one particle is equal and opposite to 

that of the other. The propogator for this system is solved exactly and expressions for 
-—, the mean displacement and mean square displacement (MSD) are found. Depending 

d , on the nature of the intrinsic bias, the system's behaviour displays two regimes, 

characterised by (i) the particles moving towards each other and (ii) away from each 
other, both qualitatively different from the case of no bias. The continuous-space limit 
of the propogator is found and is shown to solve a Fokker-Planck equation for two 
^ ' biased, mutually exclusive Brownian particles with equal and opposite drift velocity. 
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1. Introduction and motivation 

Systems of randomly moving agents that exclude one another from the space they occupy 
are ubiquitous in science and technology, from RNA transcription [H[2], to territorial 
behaviour in the animal kingdom [3] to wireless networking [4] . The theory of exclusion 
processes has been studied since 1965, when Harris [5] showed that a tagged particle, or 
tracer, on a ID line subdiffuses at long times. Since then, there have been a variety of 
mathematical developments of these so-called single-file systems [BTfTO] whereby particle 
motion is overdamped and interaction is mutually exclusive. 

When the mutually excluding particles are unbiased, one often talks about systems 
undergoing symmetric exclusion [11]. On the other hand, if the particles are subject 
to a drift, one talks about asymmetric exclusion [6]. However, in all cases studied 
so far, symmetric or asymmetric, the particles undergoing exclusion exhibit identical 
behaviours. In [12] Aslangul solved exactly a particular symmetric exclusion process: 
the case of two unbiased repulsive random walkers on an infinite ID lattice. Here, we 
extend that work to the case where each walker does have an intrinsic bias, but the 
bias of one is anti-symmetric to the other. That is, the probability of the left-hand 
particle jumping right (left) at each step is < p < 1 (1 — p) and the probability of the 
right-hand particle jumping left (right) is also p (1 — p). 

The practical motivation for our study arises from the collective emergence of 
territorial patterns in animal populations [3]. Animals are called territorial if they 
each defend a region of space from possible intruders or neighbours. Since they need 
to move around to carry out their vital activities such as foraging, animals are unable 
to monitor their territory boundaries on a permanent basis. For this reason, many 
species have evolved an ability to define their territories using scent marking, thereby 
eschewing the need for continuous border patrolling. An animal marks the terrain it 
visits by depositing a recognisable olfactory cue that is considered to be 'active' by 
conspecifics for a finite amount of time. As neighbours encounter active foreign scent, 
they move away to avoid costly confrontation. 

By modelling animals as territorial random walkers [13], that is random walkers 
with such a scent-mediated interaction process, the terrain naturally subdivides into 
territories, demarcated by the area that contains active scent. In ID, each territory is a 
finite interval joined to adjacent territories at what we call the borders. Since the scent 
is only active for a finite time, unless the animal re-scents its borders within this time, 
the borders will move. Thus the borders can be viewed as randomly moving particles 
in their own right. In addition, since smaller-than-average territories in the model end 
up having their borders re-scented more frequently than larger ones, they will tend to 
grow, whereas larger-than-average territories tend to shrink, meaning that the borders 
can be thought of as randomly moving particles connected by springs. 

In figure [IJ we sketch a mathematical representation of the territories. Each 
spring represents a territory, whose width fluctuates around a mean length equal to the 
inverse of the animal population density. Each border is a particle whose movement 
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Figure 1. Diagram of a model of territorial dynamics that reduces the interacting 
particle model of |3]. Territories are modelled as springs, as in [T3], joined together 
by borders that are modelled as diffusive particles. Zooming in on a border point 
reveals that it consists of two boundaries, each of which is moving randomly but with 
a drifting tendency towards the other. 



is intrinsically random, though also constrained by the presence of the connected 
springs. Consequently, since this is a form of symmetric exclusion process, the resultant 
movement of a tagged border particle is subdiffusive [H]. 

However, by zooming in on a border one realises that it is actually made of two 
boundaries, one for each of the two adjacent territories. The process by which the 
movement of these two boundaries gives rise to the intrinsic random movement of the 
border can be described by our present analysis of anti-symmetric random walkers and 
provides the main motivation for this work. 

The paper is organised as follows. The model description and its exact solution 
form section |2j In section |3] long-time dependences are studied and compared with 
stochastic simulations, whereas the spatial continuum limit is analysed in section HJ 
Section [5] explains in more detail the connection of this walk to systems of territorial 
random walkers and section [6] contains some concluding remarks. 

2. The model 



The starting point of our investigation consists of writing a master equation for the 
joint occupation probability P n , m (t) of the two particles being at site n and m at time t. 
In relation to the territoriality problem, the two particles are the two boundaries that 
constitute a border, disregarding the presence of other territories. Both cannot occupy 
the same site at the same time, but unless impeded by this constraint, at each hop the 
left-hand (right-hand) particle moves right (left) with probability p and left (right) with 
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probability 1 — p. The hopping rate, i.e. hopping probability per unit time, is denoted 
by F and the lattice spacing by a. As particles may only hop to nearest-neighbour sites, 
we follow Aslangul's construction [12] and write 

dP 

-^(t) = 2F[P n+1 , m (t) + P„, m _!(t)](l -p)(l - 5 n , m _!)(l - 5 n , m ) 

+ 2F[P n _ 1 , m {t) + P n>m+1 {t)]p(l - < m+ i)(l - 5 n , m ) 

- 4P(1 - 5 n , m+1 )(l - < m )(l - p)P„, m (t) 

- 4F(1 - 5 nim _!)(l - <5 n , m )pP n , m (». (1) 

The term 1 — 5 njm , where 5 is the Kronecker delta, represents the fact that two particles 
cannot hop from the same lattice site, whereas 1 — # n , m ±i represent the situations where 
both particles occupy adjacent lattice sites and so neither can move towards the other 
on the next hop. 

To seek the exact solution of ([TJ, it is convenient to use the generating function [TB] 
for P n , m (t), which is f{<j>,tj),t) = En^-ooPnAt)^^- The master equation © 
implies the following relation for the generating function 

%) = -if If ^ ^ t) - r ^ [i + COSl _ p (0 - 0')] f(4>\ 4> + 1>-4t,t) 



o^ 2tt 

' f ir /(0/ ' + ^ " 0/ ' t} [coSp(0,) + cosi -^ 



dt w ^ T,r " J 2tt 

+2P(cos p + cosi„ p ip)f((j), ip, t) 

-2F [ ^/(0', <f> + ip - 0', t) [cos p + CO Sl _ p (^)] 

n2n d0' 

'0 

where we have introduced the notation cos p (6>) = pe l8 + (1 — p)e ld so that cosi = cos. 

At time t = the particles occupy two lattice sites, denoted by A r 1 and N 2 . Without 
loss of generality, assume N\ < N 2 and since the particles cannot cross, the particle 
starting at Aq is referred to as the left-hand particle, the other is the right-hand particle. 
By using this initial condition and setting 9 = + -0, the Laplace transform of (J2J) is 

/(0, 9 - 0, e) = 0(0, 0) + V 0,(0, 0) / * A<f>%{<f>')f{<f>\ 9 - 0', e), (3) 

where /(0, 9 — <p,e) = f d£/(0, -0, t)e~ et is the Laplace transform with variable e, and 



9[ - '^ e + 2P[2-cos p 0-cosi_p(0-0)]' 

/a x\ _ 2F t 2 ~~ cos p ^ ~ cos^ p (9 - 0)] _ 1 

° l(M) " e + 2P[2-cos p 0-cos 1 _ p (0-0)]' M0) " 2V 

, fl ,, 2P(l-p)[2e^-(l+e ie )] 



e H- 2P[2 - cos p - cosi_ p (0 - 0)] ' v y 2tt 
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2Fp[2e~^ - (1 + e- w )} 



a 3 (9,<j)) 



e + 2F[2 — cosp — cosi_ p ( 



W) = 5-, (4) 



where AiV = N 2 — Ni. If, for a given value of 8, we set h(<f)) = f(<fr, 9 — 0, e) and write 
ai(4>) = ai(9, 0), g((f>) = g{6, 0) to ease notation, ([3]) can be written in terms of the single 
variable d> as follows 



h{(j>) = g{(j>) + £>(0) / WWW 



i=l 



2- 



(5) 



This is a Fredholm integral equation with degenerate kernel [15]. After some lengthy 
algebra (see Appendix A), the following solution is eventually found 



f(4>,i>,e) = e ^(^)L^ cos t±A( e u_ e (i-AN)u 



AN 



+e (l-AV)« e i^ e i(A7V-l)^±£ ( V 



.1 — p 

e + 2F(2 — cosp — cosi_ p ip) 
Here, u is defined by the equation 



_ e iA7V</> / 



1 — p 



V 



e cos 



VI — p 

<p + ?p 



i ^ -1 

P \ 2" 



1 — p 



(6) 



Z + 1 + , /(Z + l) 2 - 4p(l - p) cos 2 f 

e u = ¥ - (7) 

2[p(l-p)]?|cos|| 

where Z = e/AF and the branch of the square root function used here, and elsewhere 
throughout the text, is the one that takes real positive values when the argument is a 
positive real number. 



3. Asymptotic analysis 

In order to examine the asymptotics of the system, it is convenient to choose the initial 
conditions N± = 0, N2 = 1, as this gives rise to a simpler form for ([6]) 



f(<l>^,e) 



e ^ 1 cos 4+2 1 



■ A-iji 



4>H>\ 



4F[Z + 1 — h(cos p + cosi_ 



R{<j> + ^, Z)e* V - 2(1 - p) I cos ^ 
i?(0 + i),Z)- 2(1 - p) cos 2 S^ 



(8) 



where i?(6*, Z) = Z + 1 — -1 /(Z + l) 2 — 4p(l — p) cos 2 |. This readily reduces to a result 
of Aslangul (equation 2.11 in [12]) when p = |. 

The marginal distribution for the left-hand (resp. right-hand) particle can be 
calculated by setting ip = (resp. = 0). For — it < < it, we have the following 
expression for the generating function of the distribution of the left-hand particle in 
Laplace domain, when the right-hand particle can be anywhere else, 



/(0,O,e) 



cos I 



4F(Z + \ 



2 2 cos p ' 



J R(0,Z)e i l-2(l-p)cosf 
' i?(0,Z)-2(l-p)cos 2 | 



(9) 
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This allows us to calculate the mean position (xi(e)) of the left-hand particle in Laplace 
domain, by differentiating (Q with respect to <fr, multiplying by — ai and setting = 0, 
with the result 



(xi(e)> = £ (l - - e y / e2 + 8Fe+lQF2(l-2p)A + QF(2 g 1} . (10) 



lr 



Differentiating Q twice with respect to 0, multiplying by —a 2 and again setting = 
gives the second moment of the distribution 



a 2 
~Ae ' 






' 63 



(a^(e)) = — ( 1 + — - - Ve 2 + 8Fe + 16F 2 (1 - 2p) 2 

Uf 2 (1 - 2p) + F^e 2 + 8Fe + 16F 2 (1 - 2p)A . (11) 

By using the fact that /^[(e 2 + 2be + b 2 — a 2 ) -1 / 2 ] = e~ bt I (at), where £ _1 denotes the 
inverse Laplace transform and I„(z) a modified Bessel function of order u, expressions 
( TTUl) and ( TTTT) can be inverted exactly to give the respective formulae in time domain 



( Xl (r)) = ~ ^(2p-2)T + 8^W^ j o T ds^e-^hlS^pJT^s}^ , (12) 

( x 2 ( r )) = a 2 [{2 - 2p)r + 2(1 - 2p){2 - 2p)r 2 

+2 v ^r^)| T ds r ~ S ~ 2(1 ~ 2p)(r ~ s)2 e- 4 -J 1 [8v / Rr^) g ]), (13) 

where r = IF is dimensionless time. Denote by X\{j) and Xi{f) the positions of 
the left- and right-hand particle respectively and let d{r) = (x 2 (r) — Xi(r)) be the 
mean separation distance. Since the second moments of the particles coincide and 
(xi(r)) = — (x2(t)), it is convenient to denote by (x 2 (r)) the second moment of either 
particle and by Ax 2 (t) = (x 2 (r) — (x(r)) 2 ) the mean-square displacement. 

If p — ~ then the integrals in ([12]) and ({TBI can be computed exactly [12]. For p ^ \, 
the integrals J °° (\.ss n e~ is Ii[8^p{l —p)s\ for n = —1, 0, 1 are the Laplace transforms of 
£ n ii[8y/p(l — p)t] evaluated at the point where the Laplace variable is equal to 4, that 
is 



£{r l h[&^(l=p)t]}{e)\ 
£{/ 1 [8 v /p(l-p)t]}(e)| e=4 
C{th[8^p{l-p)t]}{ 



e=4 



1- |l-2p| 




2y/p{l-p)' 

1- |1 -2p\ 




8yjp{l-p)\l- 
\/p(1-P) 


2p| 



e " 4 8|l-2p| 3 ' 



(14) 
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Each of these three terms is finite for p ^ |, so this allows us to obtain asymptotic 
expressions for ([12]) and (TT3l) yielding the following expressions for r ~^> 1: 

if 5 < p < 1, 



d(r) 



V 



2p 



a\/T 



i 

2' 



4a(l - 2p)r 
2a 2 (l-p)r 



(x 2 (r)> 



2a 2 r 
4a 2 (1 



2_2 



2p) 2 r 



Ax 2 M 



2a 2 (l-p)r 



2a 2 (l 



2a 2 r 



7T 



if p 

if < p < |. 

if | < p < 1, 
ifp=§, 
if <p < \. 

if \ < p < 1, 

ifp=|, 
if < p < i. 



(15) 



(16) 



(17) 



The different qualitative behaviours in both the MSD and the mean separation distance 
are now evident. The limits p — > ~ and t — > oo do not commute, so the asymptotic 
diffusion constant is very different in the case P = \ from the cases where p is either 
just above or just below \. Figure |5] shows the timescales in which the three regimes 
diverge from one another. 
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Figure 2. Panel (a) shows the MSD as it varies through time for values of p close 



to 



demonstrating when the MSD begins to split into three regimes, p < i, p 



2' f 2' 

0.45, 0.49, 0.499, 0.5, 



p > ^ . Values of p from the top curve to the bottom are p 

0.501, 0.505, 0.51. Panel (b) shows the timescale t' beyond which the MSD curves for 

different values of p diverge by more than 1% from the curve for p = -i. 



For d(r), the different qualitative dependencies occur in the exponent of time so 
that for p <\ the displacement saturates, whereas for r > \ it increases. Furthermore, 
this increase is linear for p > ^ but sublinear when p = \- Figure |3] compares the various 
asymptotic expressions with simulation output for various p. 
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Figure 3. Comparison of the asymptotic expressions from (fl~5j) . (fT6|) and (fl7|) with 
average values of 10 6 stochastic simulations of the system for p = 0.4,0.5,0.6. Panel 
(a) demonstrates how the mean distance between particles d(r) exhibits qualitatively 
different behaviour in the three regions p < \ , p = \ and P > \ when plotted against 
dimcnsionless time r. Panel (b) shows the quadratic nature of the asymptotic second 
moment of a tagged particle when p > ^ , as compared with p = ^ or p < ^ when the 
second moments are asymptotically linear. In panel (c), we see the particles reaching 
their asymptotic diffusion constants. 

Conversely, at short times the behaviour of the system depends continuously on p. 
For r<l, considering only terms that are linear in r we find: 



d(r) « 1 +4o(l-p)r, 
(x 2 (r)> ^2a 2 (l-p)r. 



(18) 
(19) 



The second moment expression at short times differs from the corresponding long time 
expression by a constant for p > | but by order r for p < |. Consequently, the shape 
of the second moment's evolution over time is very different for the two regions p < \ 
and p> \i despite their identical short-time approximations (see figured]). 

4. The continuum limit 



The transition to continuous space is made by taking the limits as a — > 0, F — » oo, 
Ni — > oo, N 2 — > oo and p — )■ | such that D = a 2 F, x\$ = aNi, x 2 ,o = aN 2 and 
v = 2aF(2p — 1). Here, D represents the diffusion constant, x\ t o and X2,o the start 
positions of the left- and right-hand particles respectively and v the velocity of one 
particle towards the other, the latter of which may be positive, zero or negative. Also 
denote by Axo = £2,0 — x bo the distance between the two starting positions. 

By setting <fi = k±a and ip = k 2 a, the aforementioned limit, is found for ([6]) and 
denoted by Q(ki, k 2 , e): 
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Figure 4. Comparison of exact analytic expressions for the second moment (TH 
with short-time ([TO]) and long-time ([To 7 )) approximations. Panel (a) shows cases where 
p > i and both approximate expressions are parallel. As p increases towards 1, the 
distance between the two approximations decreases and the curves converge faster 
towards the long-time expression. Panel (b) shows cases where p < \ . The short-time 
approximations are linear whereas the long-time ones are quadratic. 



Q{h,k 2 ,e) 



e iAx k 2e ix 1 , (k 1 +k 2 ) 



e + i{k 2 - h)v + f [(h + k 2 ) 2 + (fca - fcO 2 ] 
e + i(fc 2 - h)v + f [(h + k 2 f + (k 2 - k x y] 



+ 



cxp 



Axp / u 



2D \ V2D 



^ + §(h + k 2 y 



£ + W + T(h + k 2 y- -/§ 



(20) 



2D 



This reduces to a result of Aslangul (equation 3.1 in [12]) by setting v — 0, Xi i0 = and 
£2,0 = 0. By using the identity 



C 



-1 



e -A( v / i+C- J B) 



v^+C-fi 



. P AB -4r 



7rt 



1 + erf 



2Bt- A 

2y/t 



R 2 t 



(21) 



from [IT], where erf (2) is the error function, (1201) can be Laplace inverted to give the 
following expression 



A(k 2 -k 1 )vt-§[(k 1 +k 2 ) 2 +(k 2 -k 1 ) 2 ]t I iAx k 2 



Q(h,k 2 ,t) = e - 1 ^-fciM--K*i+*2) +l*a-*i; jtJ e i^o*2 + /jgiAxo-^-*^ _ fcl ) x 



'£> 



fcl+fco 



d,s 



e sds y /Ax — zvs 

— == 1 ^=erlc - — 

7TS v^D V v / 8Ds 



-(i(fc 2 -fci)f-f (fc2-fci) 2 )s l„ ia: i,o(fci+fc2) 



, (22) 
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where erfc(z) is the complementary error function, erfc(z) = 1 — erf(z). In order to 
Fourier invert (122]) it is convenient to perform the double integral in the coordinates 
K = k\ + k 2 and k = k 2 — k\. This procedure yields the joint probability distribution 
in continuous space and time 



{xi-x 10 -vt) 2 (3=2- z 2,0+" i ) 2 



Ol-^l, 0+^2-^2,0) 



Q(x 1 ,x 2 ,t) 



VAirDt VAnDt 



VSirDt 



f l [x 2 - Xi + 2v(t - s)}e 

JO S ADy/7T(t-s) 3 



[x 2 — x 1 +2v(t— s)] 2 



SD(t-s) 



(Ax -2vs) z 
e 8D~s 



+ 



v „ /Axq — 2vs 
=ertc 



'2D 



V8D~t 



,(23) 



where Q(x\, x 2 , t) is the inverse Fourier transform of Q(k\, k 2 , t). The first summand in 
( [23l displays the short-time behaviour whereby the probability distribution of the left 
(right) particle can be approximated as a narrow Gaussian travelling right (left) at speed 
v and the interaction between the two particles is minimal. This interaction, represented 
by the second summand in (j23p . becomes more pronounced as time increases. 

It turns out (Appendix B) that (1231) is a solution to the following Fokker-Planck 
equation that is obtained by taking the continuum limit of the discrete-space master 
equation ([2]) in the region \n — m\ > 1 



dQ 



(x 1 ,x 2 ,t) = D 



2 



Q{x 1 ,x 2 ,t) + v 







Q{x u x 2 ,t). (24) 



dt ^ 1, ^ z, ' \px\ dx\) ^ x -^" 1, ^" z, " yg X i dx 2 

However, this continuum limit is only valid for X\ ^ x 2 . Since the particles cannot cross, 
and therefore the probability density along x\ = x 2 must be zero, one can interpret this 
physically by imposing a zero- flux boundary condition along the line X\ = x 2 [18], that 



is 



D 











Q(xi, x 2 ,t) + 2vQ(xi,x 2 , t) 



0. 



(25) 



x 1 =x 2 



dx 2 dx\ 

which is automatically satisfied by ( 123]) . As such, the solution reduces to a result of 
Ambjornsson et al. [18] in the case v = 0, as well as Aslangul [12] when additionally 
Xi t o = and x 2 ,o = 0. 

To find expressions for the mean separation and MSD, an identical procedure to 
the discrete case is pursued (Appendix C), giving the following results 



d(t) = Axq — vterfc 
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vt r (2vt-Ax 

— erfc ; — 

.2 V VSDl 
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as ; e « Ds 

A^2lTs 



(27) 
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where Ax 2 (t) is the MSD of either particle (Ax 2 (£) = Ax?>(t)). In the case v 
integrals in (126]) and (12TJ) can be calculated exactly to give the following 



0, the 



d(t) 



Ax 2 {t) 



— e 80?^ + Ax - Ax erfc ( Ax W^- j 
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(28) 



+ 



(29) 



For u 7^ on the other hand, the infinite integrals J„ dss n > e sd3 for n = — 1, 0, 1 
are finite, so calculating them allows us to obtain asymptotic expressions for (126]) and 
( |27l) yielding the following expressions for t ~^> 1: 



d(t) 






Ax 2 (t) 



5Di 

7T 

2ui 



if v > 0, 



if v = 0, 



if v < 0. 



(30) 



2.D(1 -~)t 

7T 

2Dt 



if 1; > 0, 
if v = 0, 
if u < 0. 



(31) 



whereby d(t) ~ Aa; — 2vt and 



This contrasts with the small-time limit t <C 1, 
Ax 2 (t) ps 2L>£ for any u. 

Notice that the ?; > (v = 0, v < 0) cases of (J3~0l and 1 J3T1) are simply the 
continuous-space limits of the p > » (p — ~, p < =) cases in the discrete-space 
expressions f fT5|) and ffTTj) . For example, in the case p > \ from (fT5l) . by setting 
a 2 r = D£ and t> = 2aF{2p — 1), we obtain d(t) = 2Dp/v and by taking the limit 
p — > I one recovers the continuous asymptotic result d(t) ~ D/v reported in ( 130]) . 
Likewise, setting a 2 r = -Dt in the case p > \ from (ITTj) and by taking the limit p — > ■= 
one recovers the continuous asymptotic result Ax 2 (t) ~ Dt from (1311) . 



5. Connection to territorial random walkers 



In [13], simulation analysis of the many-bodied, non-Markovian system of territorial 
random walkers demonstrated that the asymptotic generalised (because of single-file 
phenomena) diffusion constant of a territory border depends on an interplay between 
the so-called active scent time Tas, the time for which a scent mark is recognised by 
conspecifics as an active territory cue, and the animal population density p. Specifically, 
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the border diffusion constant decays exponentially as the dimensionless parameter 
Z = TxQ,Rp 2 a 2 is increased, where R is the rate of the animal's movement between 
lattice sites, separated by distance a. Part of the purpose of the present study is to gain 
a deeper insight into why this phenomenon is observed. 
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Figure 5. The relationship between the value of p measured from simulations of a 
system of ID territorial random walkers and the dimensionless quantity Z defined in 
secton[5j This is compared with the probability Pp(Z) that the animal fails to traverse 
a territory of average width 1/p within a time Tas- In order to measure p from the 
simulations, the number of times a boundary moved towards the adjacent boundary 
were counted, and divided by the total number of times that the boundary moved. 
The equations for the curves are 1 — p = 0.49e~ 27Z and Pf(Z) = e _7r z ^ 4 , where 
7r 2 /4 ps 2.5. The inset shows how the variance V(S) of the territory size S decays as 
Z increases, used in the main text to explain the discrepancy between the rate of the 
exponential decays of the two curves in the main plot. 



In the territorial random walk system, p is the probability that, if there is a gap 
between two adjacent boundaries, that gap will decrease in length the next time a 
boundary moves. Such a probability is clearly always greater or equal to | on average, 
otherwise the territories would fail to maintain a positive average width. For such 
values of p, equation (TTT1) shows that the asymptotic diffusion constant of a boundary 
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is proportional to 1 — p. 

When we measure the value of 1 — p directly from the simulations, we see that it 
also decays exponentially as Z increases, suggesting that calculating p is of fundamental 
importance in understanding why the border diffusion constant decays exponentially as 
Z increases. This relationship between p and Z can be explained as follows. First we 
observe that the probability of a boundary decaying is likely to be closely related to 
the probability of an animal traversing its territory within a time Tas- To this end, 
we calculate the first-passage probability ^(t) for an animal to traverse a territory of 
average length 1/p, which corresponds in our lattice system to an integer N = 1/ap 
sites. Since this is equivalent to the situation where the animal starts at a reflecting 
boundary has to traverse to the other, absorbing boundary, the asymptotic value of 
this first-passage probability is calculated in [19] to be F(t) ~ e~ n Rt ' 4:N . Therefore 
the probability Pf{Z) oc j T dtJ-(t) of failing to traverse the territory within a time 
Tas is approximately e _7r z ^. In figure El Pf(Z) is plotted alongside the simulation 
measurements for 1 — p showing that both decay exponentially with increasing Z and 
with similar exponents. 

To explain the small discrepancy in the two exponents, we make the observation 
that as Z is increased, the variance in the territory width decreases in an approximately 
exponential fashion (inset figure 4). Because of the N 2 dependence of the mean first 
passage time to cross the territory [19], the mean first passage time increases as the 
variance in the territory width increases. Therefore for a fixed p, the actual mean first 
passage time to traverse a territory decreases as Z increases, whereas above we have 
assumed that the first passage probability is always equal to that of a territory of average 
width. This has the effect of causing the probability 1 — p to decrease with Z slightly 
faster than in the analytic estimation. In other words the curve of Pp(Z) decays slightly 
slower than the curve of simulation measurements of 1 — p. 

6. Conclusions 

The propogator for a system of two anti-symmetric, biased random walkers on an infinite 
lattice is computed exactly. We characterize the bias via the parameter < p < 1, 
representing the probability for the walkers of moving away from each other, towards 
each other or with no bias at all. These three distinct physical scenarios depend, 
respectively, on the value of p being less than, greater than or equal to |. When p 
is less than |, the walkers drift away from one another with a mean displacement that 
is asymptotically linear in time. At p = | the random walkers still drift apart, although 
the mean displacement scales as the square root of time. For p > |, the distance between 
the particles saturates. The asymptotic saturation distance comes about because of two 
opposing tendencies in the walkers: a drift towards one another, given by the amount 
of bias in the walkers' movement, and the magnitude of their intrinsic diffusion. 

The corresponding propagator for the continuum limit is also computed exactly by 
setting D = a 2 F and v = 2aF(2p — 1) in the limits a — > 0, F — )■ oo and p — > \. It is 
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the solution of a Fokker-Planck equation with zero-flux boundary conditions whenever 
the particles meet. 

A motivation for this study is part of a programme to understand systems of 
territorial random walkers [3]. A first step in that direction has been the study [13] of 
a reduced one-body dynamics for the movement of a single animal within subdiffusing 
territorial borders attached by springs. The present work represents an additional step 
in the formulation of a simplified model of the non-Markovian dynamics of territorial 
random walkers, whereby the fine scale dynamics of a border are studied through the 
analysis of its constituent adjacent boundaries modelled as anti-symmetric random 
walkers with a bias probability p > 1/2. Future work will involve studying the effects 
of having a sequence of interlinked, randomly moving borders, as sketched in figure [U 
with the left and right boundary of each territory connected by springs. 
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Appendix A 

Since (J5J) from the main text is a Fredholm equation with degenerate kernel [15] . its 
solution is a linear combination of the quantities % = J^ d(p'bi((j)')h((j)') for i = 1,2,3, 
which satisfy the following system of equations 

7l(l - «ll) -72«i2 -73«13 = 01, 
-7l"21 + 72(1 - "22) - 73«23 = 02, 
~7l«31 -72«32 +73(1 +«33) = 03- (A.l) 

The various /3j's and a^'s can be calculated as t = J Q W d^'b^^g^') and a^ = 



"271 



J Q n d(j)'bi((f)')aj((f)') to yield the following expressions 
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e iiVi0 e i^ e -AiV« 



8F(1 — p) cos I sinh u \1 — p 
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8F(1 — p) cos I sinh n \l — p 
[p(l -p)]~$ -2e-"cos| 
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2 cos 2 sinh u 
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(A.2) 



In these equations u is defined by (|7|) in the main text. Solving the system of equations 
(lA.lj) eventually gives 



7i = 72 = 

eiNieei i^m e . ANu[p/{1 _ p)] ^ 



73 



4F{(l-p)cosf - [p(l -p)]2e" u } 



-2tt 



(A.3) 



Plugging these values for ji = j^ &(j)'bi(<ft)h(4>') into (E]) in the main text gives the 
expression for the generating function of the system's probability distribution in Laplace 
domain. 



Appendix B 

In order to solve ( )24l) with boundary condition ( )25l) from the main text, it is convenient 
to convert to coordinates x s = Xi — X\ and x c = (x\ + x^jl so that x s is the separation 
distance between the particles and x c is the centroid. This allows us to write (liH|) as 



9R, v „ (I d 2 d 2 \ „. , 0i2, ,„„, 

— {x c , x s , t) = D I -— + 2— 1 i?(x c , x s , t) + 2w— (x c , x a , t), (B.l, 



where R(x c ,x s ,t) = Q(xx,X2,t). The flux vector of equation (1B.1|) is 



J 



2D——(x c , x s , t) + 2vR(x c , x s , t) 
ox.. 



(B.2) 
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so the zero- flux boundary condition mentioned in the main text is n ■ J\ Xa =o = where 
n is a unit normal to the line x s = [18]. By writing R(x c ,x s ,t) = R c (x c )R s (x s ), (IB.lj) 
becomes 



and 



oR r . . Z/ c/ ix c . . 

[x c , tj = — — — r-(x c ,tj, 



Of 



2 9x2 



£/jl s / \ „ ^ d R<, , , „ oRb , . 

_ (l „ t) = 2D (x „ i)+2 „ (1<|t)i 



with the boundary condition 



(B-3) 



(B.4) 



D— — (x s , £) + vR s (x s ,t) 
ox. 



0. 



x s =0 



(B.5) 



The solution to (IB. 30 is a Gaussian and the solution to (1B.4J) with boundary condition 
(IB. 51) can be found in e.g. [20] with the result 



-**'S \Xc, £J 



(xc-k c0 ) 



V2nDt 



(B.6) 



R c {x s ,t) = H(x s 



(x s —x 3 Q+2vt) 
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V8nDt 
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v r f x s + x s0 — 2vt\ _vx 

H ertc H — e D 

2D V v / 8Dt 



(B.7) 



where x s> o = Axo and x Ci o = (x^o + X2,o)/2 are the initial conditions, and H(x) is the 
Heaviside step function (H(x) = if x < 0, iif(x) = 1 if x > 0). The solution to (J21 
from the main text can now be written down as 



(x 1 -x lt0 +x 2 -x 2i o) (x 2 -x 2: 0-x 1 +x 10 + 2vt) 



Q(xi, x 2 , t) = H(x 2 - Xi) 



V8nDt 



V2nDt 



+ 



eao 



^(X2,0-Xlfi-Vt-X2+Xl), 



(x 2 ~x 1 +x 20 -x 10 ) 
8Dt 



V2nDt 



+ 



V , (X2-Xi+X 2 fl-Xi fl -2vt 

— ertc 7 — ■ | e D 

D V y/SDi 



(B.8) 



In order to show that (IB. 81) is equivalent to (j2"3"|) from the main text, the following 
integral is calculated 
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I(x s ,t) 
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(B.9) 



Since this is the sum of two convolutions in time, its Laplace transform can be found 
by using the identity £[f * g] = £[f]£[g], where the asterix denotes the convolution 
/ * 9 = Jo dsf(s)g(t - s). Since 



[x s +2v{t-s)} 2 

[x s + 2v(t - s)]e 8 °(*- s ) 



[x a +2v(t-s)] 2 
Q e 8D(t-s) 



4Dy/ir(t-s) 3 dx s yfifii - s ) 

the Laplace transform of I(x s ,t) can be written as 
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(B.10) 



By repeatedly using the formula ( 12~TT) from the main text, expression ( IB. 101) can be 
Laplace inverted to give 



I(x s ,t) 







dx. 



V VSDi J 



Performing the differentiation with respect to x s gives 



fB.ir 



^U y s/oDt / 



sgn(x s J 



e2D 



!Dl ail > Xs OT ^P HOt 



(B.12) 



V2^Dt 

where sgn(x) is the sign of x (sgn(x) = — 1 if x < and sgn(x) = 1 if x > 0). After 
replacing the second term of (|23|) from the main text with I(x s , t) one can show that the 
continuum limit of (Q is indeed the solution of the Fokker-Planck equation (|2"31) with 
the above mentioned zero-flux boundary conditions. 
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Since the values of d(t) and Ax 2 (t) depend only on the initial condition Axo and not 
the specific values of Xio and a^o, calculations are simplified by assuming x 1)0 = 0. 
The marginal probability distribution for the left-hand (right-hand) particle in Fourier- 
Laplace domain is found by setting k 2 = {k\ = 0) in equation ( 1201) . Focussing on the 
left-hand particle gives the following expression 



Qi(h,e) 



ie iAxo ^rfci 



D\J^D V e+ 2D + 2 k l 



e - ik t v + Dk\ e - \k x v + § k\ 



6 + 2D + 2^1 



(CI) 



'2D 



This allows us to calculate the mean position (xi(e)) of the left-hand particle in Laplace 
domain, by differentiating Q with respect to fci, multiplying by — i and setting k\ — 



M*)> = 3 



£>e 



yiD v"2~D V i "2D 



^(V^+fi-vfe) 



(C.2) 



Differentiating ( 1C.1I) twice with respect to k\, multiplying by —1 and again setting 
ki = gives the second moment of the distribution 



<*)) 



2v< 



2D 

/2 
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By using the formula (1211) from the main text, (1C.2I) and (1C.3J) can be inverted exactly 
to give the respective formulae in time domain. Performing the same calculations for the 
right-hand particle allows us to find the following expressions for the mean separation 
and MSD 
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V y/8Dl 
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Ax 2 (t) = 2Dt - Ax vt + v 2 ds(t - s)erfc 
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(C.5) 



Expressions (12 6p and (T2"T|) from the main text are obtained by applying the formula 
/„* dserf {^jf) = ted (^n) - ^ / * ds^t^e~(^) throughout (JC2D and flOoD . 
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